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Abstract 

We present the calculation of coherent radio pulses emitted by extensive air showers 
induced by ultra-high energy cosmic rays accounting for reflection on the Earth’s 
surface. Results have been obtained with a simulation program that calculates the 
contributions from shower particles after reflection at a surface plane. The prop¬ 
erties of the radiation are discussed in detail emphasizing the effects of reflection. 
The shape of the frequency spectrum is shown to be closely related to the angle 
of the observer with respect to shower axis, becoming hardest in the Cherenkov 
direction. The intensity of the flux at a fixed observation angle is shown to scale 
with the square of the primary particle energy to very good accuracy indicating 
the coherent aspect of the emission. The simulation methods of this paper provide 
the foundations for energy reconstruction of experiments looking at the Earth from 
balloons and satellites. They can also be used in dedicated studies of existing and 
future experimental proposals. 
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1 Introduction 


There is a renewed interest in using the radio technique for the detection of 
extensive air showers induced by Ultra-High Energy Cosmic Rays (UHECR). 
Dedicated experimental initiatives such as CODALEMA jT] , LOFAR [2j, AERA [3], 
LOPES [4] and Tunka-Rex [5], with instrumentation far superior to that used 
in the early days of radio detection, have stimulated an impressive progress 
in the field. As a result the possibility of using the radio technique, either 
as a complement to other techniques for the detection of extensive air show¬ 
ers, or even to design new UHECR detectors based on it, is now seriously 
considered [6j. 

The serendipitous discovery of 16 fast radio pulses from air showers extending 
up to the GHz range with the first ANITA balloon born antenna system j9j 
came as a surprise. This system, operating with a frequency band between 
200 and 1200 MHz, had been originally conceived for detecting radio emission 
from neutrino interactions in the Antarctic ice sheet. However, the detected 
emission was strongly polarized in the direction defined by the cross product of 
the geomagnetic field and the shower direction, which is characteristic of radio 
signals from air showers due to the current induced by the Earth’s magnetic 
field [TlfTO] . Detailed simulations of these pulses showed recently that coherent 
emission can extend up to and above GHz frequencies [715] in spite of sim¬ 
ple dimensional arguments suggesting the contrary p2j. This is because the 
varying refractive index of the atmosphere introduces a significant effect on 
the relative travel time of emissions originating from different locations. When 
these delays are accounted for, a Cherenkov like ring can appear in which the 
signals from a large region of the shower arrive with little time delay and thus 
add coherently [8]. This angle is mostly determined by the refractive index at 
the shower maximum, simply because that is the region with more charge and 
current and hence contributing most. 

Fourteen out of the sixteen detected pulses were reconstructed to be coming 
from distinct points on the polar ice cap and showed inverted polarity with 
respect to the other two. This was interpreted as clear evidence that they were 
detected after reflection on the polar ice cap. The reflection of the radio flashes 
introduces several new aspects to the calculation of pulses at the detector that 
have not been previously addressed. Naturally reflection implies a reduction 
of the emission due to the Fresnel coefficients. The relative time delays with 
respect to detection at ground level are also altered, since the pulses propagate 
upwards after reflection towards the top of the atmosphere along a decreasing 
refractive index profile. These effects need to be taken into consideration to 
interpret the events detected by ANITA and to evaluate the acceptance of ex¬ 
periments that rely on observing radiation induced by showers from mountain 
tops, balloon payloads [T3] such as the next ANITA flight and EVA [1TJ or 
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from satellite payloads as proposed in SWORD |15] . 

In this article we simulate and describe the properties of radio pulses emitted 
from extensive air showers after reflection off a surface. Most of the calculations 
are performed in a configuration suited for a high altitude balloon flight over 
Antarctica, however the developed methods can be applied to other reflective 
surfaces and different detector altitudes. We modified the ZHAireS code [IBj to 
calculate the radio emission from air showers after reflection on a flat surface. 
We first describe the geometry and explain the assumptions made to adapt 
the program to calculate the reflected radiation. After this we generate a 
set of simulations to investigate the signal properties as a function of off- 
axis angle, frequency, zenith angle, and energy of the primary particle, and 
stress the importance of properly accounting for the Fresnel coefficients in the 
reflection, as well as for the propagation of the pulses towards the top of the 
atmosphere. In the Appendix we use a ray tracing code and a simplified model 
for the emission (that displays the main features of the predicted radiation as 
has been justified elsewhere 0123), to confirm that the effect of curved light 
trajectories can be neglected. 


2 Geometry and scales relevant for reflected events 

The appealing aspect of observing radiation of air showers after reflection is 
that a large atmospheric volume can be monitored with a single detector. 
Therefore, the most interesting geometry is given by air showers that impact 
Earth’s surface at large zenith angles. The basic geometry of the problem is 
sketched in Fig. [lj We define a rectangular coordinate system with the z-axis 
pointing upwards in the vertical direction and the x — y plane tangent to the 
Earth’s surface. The reflective surface will be approximated by this plane. The 
origin of the coordinate system is the point at which the shower axis intercepts 
the Earth’s surface which is assumed to be at a ground altitude z g above sea 
level. The zenith angle of the shower, 6 , is defined with respect to the axis. 
We define the off-axis angle 0 in Fig. [2j to describe the angular deviation of 
the emitted radiation with respect to the shower axi^] A generic detector 
is positioned at a point with vertical altitude hd . In Fig. [l] the detector is 
displayed in a special position such that it sees the reflected radiation which 
was emitted precisely along the direction of shower axis with '0 = 0°. The alti¬ 
tude at which shower maximum (X max ) is reached, h-xmax, is also of relevance. 
Besides determining the angle at which the emission is largest [8j, it also sets 
the scale of distances the pulse has to travel to reach the detector, dx max + dd, 
where dx ma x and dd respectively denote the distances from the origin of the 

1 0 as depicted in Figgis also used later in this work to refer to the location of 
the observer. 
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Detector 



Fig. 1. Basic geometry for reflected signals from air showers (see text for details, see 
also Fig. [2]). 


Antenna 



Fig. 2. To describe a location of an antenna we use the “off-axis” angle ip, as that 
between the shower axis and the line joining the location of X max and the position 
of the antenna. 

coordinate system to shower maximum and to the detector. 

We illustrate the typical scales of this geometry by showing in Fig. [3] some 
parameters for a high altitude balloon over Antarctica. We take a detector 
at a typical altitude of h d = 35 km, a ground altitude of the ice cap of z g = 
2 km and the Earth’s radius Re = 6357 km. The distance d d becomes simply 
a function of 9 which is illustrated in the top panel of Fig. [3j To estimate 
the distance to shower maximum, dx max, we use the average slant depth of 
shower maximum (X max ) as observed by the Pierre Auger Observatory [I8lfl9j 
together with the atmospheric density profile used in the air shower simulation 
package AIRES [20]. Clearly the average position of X max depends on shower 
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energy, but the effect is accidentally small according to the measurements 89 
at the Pierre Auger Observatory which indicate little change of A max in the so 
energy range of 10 17 ' 8 to 10 19 - 6 USES]. The results for different primary energies 91 
are shown as a function of 9 in the middle panel of Fig. [3j It should be 92 
noted that the measured RMS fluctuations of A mav are between 20 and 60 93 

g cm -2 jl8UT9] corresponding to variations of dx max below 11 % (3%) for a 94 
zenith angle of 60° (85°) degrees. The variation relative to the total distance 95 


travelled by the pulse dx ma . Y + dd reduces to 1.5% (1.0%). 96 

A relevant parameter is the Cherenkov angle, i/)c ~ cos _ 1 (l/n), at the location 97 
of X" max that is directly obtained using the refractive index n at the correspond- 98 
ing altitude, hxmax • In the following the refractive index is approximated by 99 
a simple exponential function of altitude h given by n(h) = 1 + r] 0 e~ Kh , with 100 
r/o = 325 x ICR 6 and k = 0.1218 km -1 [15]. The Cherenkov angle at X max is 101 
shown in Fig. [3] (bottom) for the above parameterization. 102 

3 Introducing reflection in ZHAireS 103 


ZHAireS [[16] is a simulation program that has been developed combining 104 
the AIRES package for air shower simulations [20] with the “ZHS algorithm” 105 
which was originally developed to calculate radio emission from high energy ioe 
showers in homogeneous ice [2TT122] and then extended for use in air show- 107 
ers [I61I23] . The contribution from each track element to the radio pulse at ios 
any given position and time is calculated in ZHAireS assuming the particle 109 
travels at a constant speed and in a rectilinear motion and also accounting for no 
the travel times taken by the radiation to reach from each of the ends of the m 
track mm- To calculate travel times in ZHAireS we perform a numerical 112 
integration to account for the variation of the index of refraction with altitude 113 
assuming that the emission travels in straight lines to the observer (see below), m 
The curvature of the Earth’s atmosphere is fully accounted for in AIRES. 115 

We have modified the ZHAireS code to deal with the reflection of air shower ns 
radio emission on a surface. For each rectilinear track element we first find m 
the point on the reflection surface and the angle of emission of radiation with ns 
respect to the track so that the emitted ray propagates first to the reflection 119 
point and then upwards towards the observer at a fixed position. Approximat- 120 
ing the reflection surface to a plane makes it trivial to obtain this reflection 121 
point for a ray coming from any point in the atmosphere. Once this is known, 122 
the time delay due to the refractive index is easily calculated integrating the 123 
travel time over the total path of the ray before and after the reflection. We 124 
assume that the emission travels in a straight lines to the observer. We have 125 
explored the validity of this approximation using a simple one dimensional 126 
model to produce pulses which resemble the fully simulated ones. This model 127 
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z„ = 2 km, h = 35 km 

& d 





Fig. 3. Top panel: distance from the reflection point to the detector (d^ in Fig. [TJ; 
rniddel panel: distance to shower maximum (dxmax in Fig. [I]); and bottom panel: 
the Cherenkov angle at X max as a function of shower zenith angle. For the location 
of X max we used the average value of (A max ) as a function of energy as measured 
by the Pierre Auger Collaboration mm- 
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allows the calculation of travel times using numerical methods to propagate 
signals along curved trajectories. The properties of the pulses obtained with 
straight ray approximation are equivalent to those accounting for curvature 
up to shower zenith angles 9 ~ 85° and frequencies of ~ 1 GHz. This study is 
described in detail in the Appendix. 

As mentioned before, we approximate the reflection surface to the xy-plane 
defined in Section [2j assumed to be perfectly flat. The bulk of the emission 
has been shown to be concentrated in a cone that makes a small “off-axis” 
angle -0 to the shower direction as shown in Fig. [2} This angle is very close 
to the Cherenkov angle (0.5° — 1°) at an altitude at which shower maximum 
occurs |8j. As a consequence the illuminated region on the reflective surface is 
relatively small, of order 0.5 km x 1 km (1.5 km x 10 km) for a 9 = 60° (80°) 
shower. As a result it is reasonable to ignore the differences in the orientation 
angle and altitude of the reflecting surface at the locations of the different 
reflection points across the illuminated area. These differences are below 0.1° 
and a few meters respectively even for showers of 9 — 80°. There are other 
important aspects to the flat mirror approximation. When rays are reflected on 
a convex and rough surface they will diverge after reflection, and therefore the 
power received at a given surface element will be typically less than when a flat 
reflector is assumed. These effects have been studied in [15] and [25], and have 
little impact at moderate zenith angles, however they can become significant 
for high zenith angles (6 > 80°). Results of the simulations with the reflective 
surface assumed flat can be corrected “a posteriori” following the procedures 
outlined in [T51I25] . Such calculation is however very detector specific, and 
out of the scope of this article. Despite this, we expect the simulation method 
presented in this work to be very suitable for that purpose. 

At each reflection point the Fresnel coefficients are applied to the time-domain 
electric field to calculate the attenuation of the components with polarization 
parallel, ry, and perpendicular, r±, to the reflection plane, defined by the 
normal to the reflecting surface and the direction of the radiation: 


r _l 


ri\ cos 9 — n 2 ^ 


sind) 2 

rii cos 9 + n 2 ^ 

A- 

sind) 2 


and ry 



1 

(— sinfA 

\ n 2 J 
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A- 

| — n 2 cos 9 



( ^ sin 6' 

\ n 2 t 

2 

ni\ 

A- 

) + ri ‘2 cos 9 


( 1 ) 


The Fresnel coefficients for an air-ice interface are shown in Fig.|4]as a function 
of zenith angle 9 of the incident ray. A large fraction of the components of 
the field is not reflected below 9 ~ 70° — 80°, and in fact at the Brewster 
angle at 9 ~ 53° the parallel component is not reflected at all. The coefficients 
change rapidly above 9 ~ 60°. Clearly they have a drastic impact on the the 
overall amplitude, the polarization and the zenith angle dependence of the 
radio signal as will be shown in Section [4] 
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Fig. 4. The Fresnel coefficients for an air-ice interface with refractive indices 1.0003 
and 1.31 respectively, as a function of the zenith angle 6 of the incident ray. 

4 Simulations for a high altitude detector 

4-1 Simulation set and considerations 


In this section we apply the implementation of the reflection in ZHAireS to a 
configuration that resembles a high altitude balloon experiment over Antarc¬ 
tica, such as the ANITA flights, or the planned EVA mission m 

We place antennas at a fixed altitude ha = 36 km above sea-level, and choose 
the reflecting surface to be at z g = 2 km above sea-level. We adopt a refractive 
index of n = 1.31 [26] consistent with ANITA measurements of the reflected 
image of the Sun [27]. The geomagnetic field is chosen to have a typical valucf 2 ] 
of 55 /iT and an inclination of —70°. We generated proton showers with zenith 
angles 9 = {57°, 64°, 71°, 78°, 85°} and azimuth such that they always arrive 
from the geomagnetic west. For each zenith angle, we generate air showers 
with energies E = {10 178 ,10 184 ,10 19 ,10 196 } eV. We select simulations that 
have a X max similar to the average (AAax)Auger observed at the Pierre Auger 
Observatory. To do so, we pre-simulate seven air showers per configuration 
with different random seeds and we select the air shower closest to (A max )Auger- 
This results in an average deviation of |A max —(A max ) Auger | ~ 18gem' 2 , which 
is within the root mean square of the energy-dependent V max -distributions 
that have been observed. The shower simulation is run with AIRES using 
QGSJETII.03 hadronic model interactions with a thinning level of 10 -5 [S0J. 


2 See for example http://www.ngdc.noaa.gov/geomag/geomag.shtml 
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Fig. 5. Top left panel: distribution of the radio signal (flux density <b) as a function 
of off-axis angle 0 (Fig. [2]) and frequency / for an air shower with 9 = 71° and 
log w (E / eV) = 18.4. In the bottom left panel we show the distribution of the radio 
signal as a function of off-axis angle at various frequencies. In the right panels we 
show the radio signal distribution as a function of frequency, in the top right panel 
for off-axis angles equal or smaller than the Cherenkov-angle at the X max of the 
shower i.e. 0 < 0.7°, while in the bottom right panel for 0 > 0.7°. 


4-2 Results 


181 


To illustrate some of the typical features of the radio signal, we display in Fig. [5] i 82 
the flux density $ as a function of frequency and off-axis angle 0 for an air i83 
shower with zenith angle 9 = 71° and a energy E = 10 17 ' 8 eV. The flux density m 
in this case is defined as the power spectrum at a fixed frequency / averaged iss 
over a period of 10 ns, and is given in units of pW rn~ 2 MHz -1 throughout ise 
this paper. In the top left panel of Fig. [5] we show the two-dimensional dis- i87 
tribution as a function of 0 and / which displays coherent properties and is iss 
clearly beamed around the Cherenkov angle at ~ 0.77°. This can be better i89 
appreciated in the bottom left panel where we show the off-axis distributions wo 
for different frequency components of the pulse. As the frequency increases m 
the radiation adds coherently only within a smaller angle off the Cherenkov 192 
cone. In the right panels in Fig. [5] we show the spectral shape of the flux 193 
density for a variety of observation off-axis angles. At very low frequencies 194 
(/ < 10 MHz) the flux density increases until it reaches a maximum in the 195 
range (/ ~ 10 — 150 MHz) and then decreases with an exponential fall-off we 
to first order. A very important feature is illustrated in the right panels, the 197 
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Fig. 6. Flux density $ at a frequency / = 300 MHz as a function of off-axis angle if 
as obtained extrapolating the ZHAireS simulated signal at ground to the detector 
(black dots and solid line), and simulating the reflection as explained in the text 
(red crosses and dashed line). The flux density is shown for showers of E = 10 18 ' 4 
eV and two zenith angles 0 = 57° (left panel) and 6 = 71° (right panel). The 
Fresnel-reflection coefficients are accounted for in all cases. 

steepness of the fall-off has a clear dependence on the off-axis angle if of the 
detector. This dependence is key to the energy determination of UHECRs with 
ANITA as pointed out in [28]. 


4-2.1 Implications of the reflection 

Other efforts to simulate reflected radio signals from air showers have relied on 
pulses simulated at ground which were extrapolated using the attenuation of 
the signal with increased distance (\E\ oc 1/r) after accounting for the loss of 
signal induced by the Fresnel coefficients [28]. In a homogeneous medium this 
“specular approach” can be expected to be a good approximation provided 
the pulse can be considered to be in the Fraunhofer limit. As a result it can be 
expected to work better for highly inclined showers, since the distance between 
the observer and air shower increases as the zenith angle rises. 

In this work all track contributions to the pulse are reflected at the interface 
to account for attenuation with distance, for the Fresnel-reflection coefficents 
attenuating the parallel and perpendicular components of the field, and for 
the fact that reflection also alters the relative time delays of emission from 
different regions of the shower affecting the coherence properties of the pulses. 
Therefore, this method can also be applied when the reflector is not in the 
Fraunhofer limit. 

We compare the specular approximation to the results of the full ZHAireS 
simulation including reflection to illustrate the difference between the two 
methods. For this comparison we evaluate the flux density at / = 300 MHz 
at a few ground locations scaling it to account for the distance from ground 
to the location of the high altitude balloon. In Fig. [6] we display the flux 
density as a function of if for two different zenith angles. We note that for 
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Fig. 7. Comparison of the flux density $ at a frequency / = 300 MHz as a function 
of the off-axis angle i/j before (left) and after applying Fresnel reflection coefficients 
(right). Different sets of curves correspond to different shower zenith angles 8 as 
labeled. The simulated showers have an energy E = 10 18 ' 4 eV. 

6 = 57° (Fig. [6] left) the distribution in i/j is significantly wider what can lead 
to orders of magnitude of over-estimation of the flux density for the larger off- 
axis angles. For 6 = 71° (Fig. fright) we still see relevant deviations between 
the two methods, but they are significantly reduced compared to the lower 
zenith angle case. Significant deviations are also found at other frequencies. 
Moreover, the shapes of the frequency spectra obtained with the two methods 
also differ appreciably. 

It is interesting to explore how the radiation changes with zenith angle for 
a primary particle with fixed energy. As the zenith angle increases the flux 
density $ decreases. This can be seen in the left panel of Fig. [TJ. displaying $ 
at / = 300 MHz as a function of ifj for an air shower induced by a primary 
particle with energy E = 10 18 ' 4 eV. The dominant effect in the decrease is the 
increasing overall distance to the detector with 6 (see Fig. [3]). Other effects 
however compensate the decrease in <f>. The angle a of the shower axis to the 
Earth’s magnetic field at the South Pole increases in the range of 9 shown 
in Fig. [7j and the geomagnetic contribution is known to scale with sin a. 
Also showers of increasing 9 develop in a less dense atmosphere where the 
geomagnetic contribution to the electric field is expected to be increasingly 
larger [3D]. The net result, including other more subtle effects [2D], is a decrease 
of <P with 9. 

To illustrate the importance of accounting for the Fresnel reflection coefficients 
they were artificially set to 1 in the simulations shown in the left panel of 
Fig.0 while in the right panel they are taken into account. Comparing both 
panels, the peak value of the flux density is largest at relatively high zenith 
angles (9 ~ 80°) when the Fresnel coefficients are accounted for, contrary 
to what is seen in the left panel where the peak value of $ is achieved at 
the smallest zenith angles. This suggests that detection can be expected to 
be most favorable for 9 around 80°. A thorough calculation of the acceptance 
integrating over area and solid angle [T3J should also account for the reduction 
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Fig. 8. The flux density at frequency / =300 MHz as a function of the energy of 
the primary particle for three off-axis angles in a shower of 6 = 71°. The results of 
fitting a straight line log 10 = alog 10 (T7/eV) + b are shown. 

of the Cherenkov angle as the zenith angle rises (see Fig. [3]), and for the 
directionality of the detection system. Such calculation is out of the scope of 
this article. 


4-2.2 Energy dependence 

From the set of simulations we have examined the energy dependence of the 
radio signal. As before, we use the flux density at a reference frequency / = 
300 MHz for a shower of 9 = 71°. 

In Fig. [8] we select three off-axis angles and plot the flux density as a function of 
the primary particle energy. We fit a simple linear function to the dependence 
of log 10 <f> on log 10 (£) and find a slope that is consistent with 2. This confirms 
that the received flux density scales quadratically with the primary energy and 
the amplitude of the electric field scales linearly with it. This is not surprising 
since for a coherent signal it is expected that the amplitude of the electric 
field scales with the number of electrons in the shower which is proportional 
to the energy of the primary particle. We verified that the quadratic relation 
between flux density and energy of primary cosmic ray particle holds for all 
the zenith angles and all considered frequencies in our simulation set. This 
has important consequences since measuring the flux density at a given off- 
axis angle provides a measurement of the energy of the air shower. In practice 
the off-axis angle can be related to the exponential fall-off of the flux density 
as can be seen in Fig. [5j This relation is key to the energy determination 
of UHECRs with ANITA [28]. This means that it is in principle possible to 
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deduce the primary energy from a single location as long as the exponential 
fall-off of the spectrum can be determined. 


5 Conclusions 


In this paper we implemented the treatment of surface reflection for radio 
signals from air showers, upgrading existing ZHAireS simulation^" 3 ] 

As a case study we simulated radiation from a set of air showers at a high 
altitude position over Antarctica inspired by the ANITA experiment which has 
provided the only measurements of reflected radio pulses from air showers up to 
now. We described and explained the flux density distributions as a function of 
frequency, energy and zenith angle of the primary cosmic ray particle. We have 
stressed the importance of accounting for the Fresnel-reflection coefficients 
attenuating the components of the held, as well as accounting for the fact 
that reflection also alters the relative time delays of emission from different 
regions of the shower affecting the coherence properties of the pulses. 

A clear quadratic correlation between the radio observable (flux density) and 
the shower energy has been observed in ZHAireS simulations with reflection. 
The intercept of the correlation depends on the off-axis angle which can be 
related to the exponential fall-off of the spectral distribution of the flux density. 
This provides the basis for the determination of shower energy from a single 
location [281] . 

Several approximations were made in the implementation of the reflection in 
the ZHAireS code. The reflective surface is assumed to be a plane and we 
argued that this is a good approximation as long as the shower zenith angle 
is 6 < 80°. Also the emitted radiation is assumed to travel along straight 
lines to the ground and then to the detector, an approximation that has been 
extensively verified in Appendix A. 

Finally, we foresee many applications of this code to the design and feasibil¬ 
ity studies of future large exposure cosmic ray detectors based on the radio 
technique applied to air showers [T5IIT3] . 


3 Code available upon request. 
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A Appendix 354 

A.l Straight vs curved ray propagation 355 

The variation with altitude of the index of refraction of the atmosphere is 355 
known to induce curvature in the path of the radio waves when propagating. 357 
It has been shown that the effects are negligible for most shower geometries 35 s 
and observers on ground EH. for which the propagation along straight paths 359 
is a good approximation. In the case of showers at large zenith angles and 3eo 
especially when accounting for reflection, the involved distances from emission 36i 
to the detector become large (see Figs. [I] and [3]), and the curvature of the rays 362 
can be expected to increase. 353 
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To evaluate if the approximation of straight light ray propagation still holds in 
the typical geometries involved in reflection, we have developed a simple ray 
tracing code. We divide the atmosphere in many layers with constant distance 
between them. The layers are taken sufficiently narrow so that the ray can be 
approximated as traveling in a straight line along a constant refractive index n 
in each layer given by the exponential model in Section [2j The ray is refracted 
at each interface, taking into account the different refractive index at each 
layer. The travel time of the ray is calculated as the sum of the times it takes 
to cross each layer. When accounting for reflection on the ground, the ray is 
also propagated upwards through a decreasing refractive index profile until it 
reaches the detector. The arrival time assuming straight line propagation to 
ground and then to the same detector position is also calculated. The reflection 
surface is assumed to be at sea level for these calculations. Since the gradient 
of an exponential atmosphere is largest at sea level, it can be expected that 
curvature effects for reflection from surfaces at higher altitudes will have less 
impact than estimated here. 


In Fig. A.l (top panel) we show the relative arrival times of radio signals 
emitted from different positions along the shower axis of a 9 = 70° shower. 
They have been calculated with the straight and curved ray approximations 
for two particular observer positions such that the radiation arriving from 
shower maximum makes an off-axis angle ~ 0.77°. One observer is located 
on the ground and receives the rays directly, the other observer is placed at 
an altitude h ( j ~ 33 km and receives the rays after reflection on the ground 
(see Fig. [Tj) . The difference in the arrival times between curved and straight 
ray propagation is hardly noticeable in the scale of Fig. |A.1| It is in fact 
below ~ 50 ps, which corresponds to a frequency of ~ 5 GHz when using a 
quarter wavelength criterion for coherence. As a result we expect the straight 
ray approximation to be valid below this frequency. 


Similarly in Fig. A.l| (bottom panel) we show the arrival times of the pulses 
for a 9 — 85° shower and observation at an off-axis angle = 0.4°. Although 
in this case the difference between curved and straight ray propagation is 
sizable, it is an almost constant offset along shower development of ~ 0.3 
ns for the observer on the ground and ~ 0.9 ns for the observer at high 
altitude hd = 50 km. This global offsets induce unobservable phase shifts 
in the field at the detector. When accounting for these offsets, the relative 
differences between the straight and curved propagations are ~ 200 ps for the 
observer at ground level corresponding to a frequency of ~ 1.25 GHz using a 
quarter wavelength criterion, while for the observer at hd = 50 km altitude 
the differences are below ~ 20 ps (~ 12 GHz frequency). 


In the top and bottom panels of Fig. A.l the relative arrival times at ground 
and at the high altitude observer are approximately flat within a large region 
(~ 10 km) around shower maximum. As a result the emission from this region 
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70° shower 



85° shower 



Coordinate along the shower axis [km] 


Fig. A.l. Top: Relative arrival time of rays emitted along the axis of a 9 = 70° 
shower at two particular antenna locations: one located on the ground (red lines) 
and another after reflection towards a high altitude balloon (blue lines) both at an 
off-axis angle close to the Cherenkov angle t/r ~ 0.77° (see Fig. [2]). The straight 
ray approximation (dashed lines) and the curved ray propagation (solid lines) are 
shown (see text for details). A 10 19 eV, 6 = 70° Gaisser-Hillas shower profile is 
superimposed (black line). Bottom: the same as in the top panel, but now for a 
9 = 85° shower viewed at an off-axis angle ij} ~ 0.40°. 
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is coherent up to GHz frequencies. An interesting feature in the top panel 
of Fig. A.l is the inversion of the arrival times at the high altitude observer 
position. The radio signal emitted in the Cherenkov angle (in this particular 
case from the region around shower maximum) arrives last at the high altitude 
location, contrary to what could be expected in a homogeneous medium. The 
time inversion is also seen for the observer on the ground for the 85° shower 
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in Fig. A.l| (bottom). The effect does not seem to have relevant implications 
for detection. 


A. 2 A simple model to describe the pulses 


The comparison of the travel times for the straight and curved ray calculations 
already sheds light onto the validity of the straight ray approximation. How¬ 
ever a comparison between the pulses obtained with the straight and curved 
approaches in a simplified calculation will give us an estimate of the quanti¬ 
tative uncertainty in the properties of the signal at the detector when using 
the straight path approximation. 


For that purpose, we model a shower with zenith angle 9 as a one-dimensional 
charge distribution varying with time N(t) as the shower propagates along a 
given direction at the speed of light c ra. In the Fraunhofer approximation 
[23] , the electric field induced by this line of charge is proportional to the 
charge N(t), to r _1 (t), with r(t) the distance between the emission point and 
the observer position, and carries a phase factor to account for the time delay 
between the arrival of the signal emitted from different positions along the 
line: 

/ e iLot a (t) 

d t (a- 1 ) 

Here uj is the angular frequency and t a is the arrival time at the observer of a 
signal emitted at time t, where 

77 

ta(t)=t + -r(t) (A.2) 

c 

This approximation was developed for a homogeneous, isotropic and non- 
conductive medium but it can be extended to the atmosphere accounting for 
the time delays from propagation of the rays in the altitude-dependent index 
of refraction [ T6] . 


In Fig. |A.2 the pulses obtained in the full ZHAireS simulation modified for 
reflection are compared to those obtained with Eq. (A.l) using as input a one- 
dimensional Gaisser-Hillas distribution, N(t), with the same depth of shower 
maximum as the simulated shower. The amplitude of the reflected field is 
shown for observers at a high altitude located at different off-axis angles ip (see 
Fig.@> and constant overall path distance. The results have been normalized 
to the peak of the electric field as predicted with Eq. (A.l). The shape of the 
angular distribution is well described by the simple model in a wide frequency 
range. It should be noted here that, being one-dimensional, the approach can 
not fully reproduce the frequency spectrum as obtained in the full simulations 
near the Cherenkov angle, where the lateral spread is of utmost importance 
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436 


[16] . We are however confident that the approximation is sufficient to test the 
validity of the straight ray approximation. 437 




Fig. A.2. Top: amplitude of the Fourier transformed electric field as a function of 
the off-axis angle for different frequencies. Results from full ZHAireS simulations 
and the model in Eq. (A.l) are shown. The shower has E = 10 19 eV, 9 = 60° and 
a height of A max above ground hx max ~ 4.4 km. The observers are placed at a 
constant distance from X max of ~ 83.3 km. From top to bottom the observation 
frequencies are 50, 300 and 900 MHz. Bottom: Same as in the top panel, but for a 
9 = 80° shower with hxmax ~ 13.5 km and with the observers placed at a constant 
distance of ~ 276.9 km. 


A.3 Validity of the straight ray approximation 


We used the simplified one-dimensional model described above to test the 
straight ray approximation. The arrival times t a at the detector to be used in 
Eq. (A.l) are calculated integrating the travel time along both straight and 
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curved paths using the ray tracing algorithm explained above. The electric field 


is obtained numerically, discretizing the integral in Eq. (A.l) in such a way 


that the time intervals correspond to regions of the shower with approximately 
constant N(t), r(t ) and t a {t). 

The results using the straight and curved ray calculations are compared in 
Fig. |A.3 where we plot the modulus of the electric field as a function of the 
offset angle of the antenna ip for showers of 9 = 70° (top) and 9 = 85° 
(bottom). As anticipated from arguments concerning the travel times discussed 
at the beginning of this Appendix, the difference in the angular distribution of 
the electric field between the straight and curved ray propagation is negligible 
for the 6 = 70° shower at all the frequencies we probed. It can be appreciated 
that even for 9 = 85°, the effect is still negligible up to a frequency of 900 MHz. 
This justifies the straight ray approximation for the calculation of the angular 
distribution of the field, even at high zenith angles and up to GHz frequencies. 
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Fig. A.3. Top: Electric field modulus as obtained with the simple model in Eq. (A.l), 
as a function of the off-axis angle ip for three frequencies. The observers are located 
at different ip angles on the ground (red solid lines) and at an overall distance 
~ 130 km to A" max after reflection (blue dashed lines). The shower has 9 = 70° and 
I m« at an altitude above ground hx max ~ 10.4 km. The results of curved (points) 
and straight (lines) rays calculation are shown. From top to bottom, the observation 
frequencies are 50, 300 and 900 MHz. Bottom: the same as in the top panel with 
the observers located on the ground and at an overall path distance of ~ 648.5 km. 
The shower has 9 = 85° and hx max ~ 16.5 km. In both panels the fields on the 
ground are rescaled with the corresponding ratio of the distance to the ground and 
the total distance to the detector for visibility. 
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